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This docMsgnt ooBtaiins a ^^ign a^^ificaticm for softimre to 
iaplcBont a ■nlti- 1 csp oral Bay» classifier. This classifier 
uses ■ttlti-tcBporal piacel data in cosbinati<w vith transition 
probabiliti^ to correct for Misregistration and crt^ planting 
changes. 
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laplesientaticm will be on the Purdue-LARS coo^uter system. It 
will be incorporated into the ^X>-LARS¥S systoi as an add-on 
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2. APPLICABLE DOCUMENTS 


• Pinal Design Specification for ECHKLARSYS Procedure 1, Houston, 
Texas, August 1977, JSC>12742, LEC-10417. 

e Job Order 63>1347>1695 

e TIRE 


e J. R. Welch and K. G. Salter, "A Context Algorithm for Pattern 
Recognition amd Image Interpretation", IEEE Trans. SML, 

Vol. SMC-1, No. 1, pp. 24-30, January 1971. 



3. SYSTEM DESCRIPTION 


3.1 HABDWARE DESCMPTIOM 
M/A 


3.2 SOFTWARE DESCRIPTIOM 

The Multi-teiqporal Bayes classifier will proceed along the 
following cosqputational lines. 


Input Cards: Control card data, including number of acquisitions, 

assets of acquisitions for classification output 
and errar output, sun angle correction, unit and 
file nundiers, channels selected. 

Pixel Fields to be classified. 

A-priori probabilities of categories. 

Transitional probabilities. 

Test fields for error computation. 


Input Files: 


A multi-acquisition MSS data tape. 

Set of conditional probability density estimates 
f(z^^^jcj), j = category number, i = acquisition 
number, for pixels z contained in the pixel fields 
on MAPJAP or DSTTAP files. (MAPTAP input will 
provide the various inputs needed to compute 


f (z 


(i) 


Cj) as a multi-normal distribution.) 


The pixel fields to be classified must be contained in the fields 
on MAPTAP or DSTTAP. Also, the number of acquisitions available 
from MAPTAP or DSTTAP must be at least as large as the number 
called for under the card input. 


Classification will proceed according to the method described in 
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this section. Results will be outputted in EOD-LARSYS MAPTAP 
file format for selected subsets of acquisitions. 

Error computations will proceed according to the method to be 
described in this section. Printer output of these computations 
will be presented for selected subsets of acquisitions. 

Prior to classification, the estimates 

{f(z^^Mcj)) i=l, NA, j=l, NOCAT 
will be extracted from DSTTAP files or MAPTAP files. 

In the case of DSTTAP files provided, there must be at least NA 
files available. The channels specified in control card CHANNEL 
must match those specified in DSTTAP files. 

In the case of MAPTAP files, the estimates {f(z^^^ Ic^)} will be 
built from multi-normal information (covariance matrices, mean 
vectors) . The channels specified in control card CHANNEL must 
be included in the MAPTAP file aggregate. 


COMPUTATIONAL PROCEDURE FOR CLASSIFICATION: 

For classification on the first acquisition of any pixel in 

the test field described under (3), select the category j which 
maximizes 


Flj (z)=f (z^^^ 



) 


For 2 acquisitions, select the category j which maximizes 

( 1 ) 


F2j (z) = (: 


I (2) 

ic^ )f(^ 


c.(2))p(,.(2) 









For 3 acquisitions, select the category k to maximize 


P(c <3'|c '^>o ‘^')P(c <1') 

R, J X J X X 


For 4 acquisitions, select the category £ to maximize 


r4,.zi^y;r’rf(z<»ic,'^>)f(z<^>ic,<^’if,z<^>ic,<^i;» 

1 j k 




P(c '3>lc.'l'c.'2')P(c.<2'|c.'l>)P(c.<ll) 
k *x 3 3 1 1 


Output MAPTAPs are required. 

Error Calculations 

These will be based on the test pixel fields. 

Bayes risk option: 

Compute the number of test pixels N. Then, for a particular 
number of acquisitions NA, NA=1 , 2, 3 or 4 , an estimate of the 
Bayes risk is 


Rbe(NA) 


1 N 
- y- 

N 


FNA ( z . ) 
£ ^ 


FNA (z^) 


z. refers to I z. 

X I 1 


(NA) 

X 




Misclassif ied Pixel Count Option; 

Compute the number of test pixels N. For the NAth acquisition, 
compare the classification (input) of each test pixel with that 
of the NAth assignment made by the multi- temporal Bayes procedure 
described above. 

Count the pixels where this differs. Dividing this count by N 
yields Rgj,(NA). 

3.2.1 LINKAGES 

The Multi-temporal Bayes classifier will use the FORTRAN IV-G 
compiler, IBM system routines, EOD-LARSYS utility routines and 
EOD-LARSYS common blocks GLOBAL, ISOLNK, and INFORM. A new 
common block MULBAY will be provided for transfer of information 
among subprograms. 

3.2.2 INTERFACES 

The Multi-temporal Bayes classifier will require an MSS input 
file for pixel radiance value input. 

It will be invoked from the MONTOR routine of EOD-LARSYS upon 
encounterin a control card of the form $MULBAY. Back-to-back 
runs with the DISPLAY processor of EOD-LARSYS will be possible. 

3.2.3 INPUTS 

a. Processor Card 
Keyword 


$MULBAY 






b. New Control Cards 


Keyword 


Parameter 


Function 


CATNAM 


category name^^,... Names of categories. These 
category name^^Q^^^ names must match those on 

MAPTAP or DSTTAP. 


CHANNEL 


^1' * • * '^N0FET2 


Selected channels. 


NCPASS N 

(Default; N=4) 


Number of channels per pass. 


NOTE: From this channel information, the number 

of acquisitions NA will be computed. 


SUNANG 


TAPE 

or 

"l' • • * '^NA 


If tape specified, sunangles 
will be extracted from the 
header of the MSS data tape 
(Universal format) . 


NASUB 


DATAFI 


FORMAT 

MAPTAP 


MAPTAP 


ni, . . • /Hgyg 


Subset of integers 1,2,..., 
NA for print of classifica- 
tion maps and error results. 


INPUT/UN IT=N,FILE=M MSS Data Tape unit and file 
(Default: N=3,M=1) numbers. 

(Not required if 
DSTTAP input specified) 


U or L Universal or LARSYS Unit 

OUTPUT/UNIT=N,FILE=M number and beginning file 
(Default; N=2,M=2) number of classification 

map file. 

INPUT/UNIT=N,FILE=M Unit number and beginning 
(Default; N=2,M=1) file number of classification 

map file. 


Keyword 


Parameter 


Function 


DSTTAP 

INPUT/UNIT=N , FILE=M 
(Default: None) 

APRIOR 

^1' • * * '^NOCAT 
(Default: None) 

ERROR 

BAYES 


or 


COUNT 

*END* 



Unit number and beginning 
file number of classification 
map file. 

A-priori values for each 
category 

Select error computation. 


Terminates control card 
input. 


c. Field Cards 

Field cards will be handled as in EOD-LARSYS, except that 
DESIGN header cards will include SAMPLE and TEST designators. 

Sample fields for classification will be user input. Test 
fields for error computation will be user supplied (optional) . 

3.2.4 OUTPUTS 

MAPTAP files will be outputted. 

3.2.5 STORAGE REQUIREMENTS 

TBD 


3.2.6 DESCRIPTION 


The following list of subprograms will constitute the building 
blocks of the Multi-temporal Bayes classifier. 


Subprogram 

MULBAY 

SET16 

RDDAT 


EXTRAC 

CLSSAM 

CLSTST 

ERRBAY 


Function 
Main driver. 

Will read control cards. 

Will read pixel fields (sample and test) , 
establish field vertices, extract pixel 
radiance values from MSS data tape, perform 
sun angle correction. 

NOTE: If DSTTAP input, radiance values are 

not required. 

Set up {f (z |Cj ) } 

Will perform sequential classification of 
sample fields. 

Will perform sequential classification of 
test fields. 

Will compute Bayes error. 
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Subprogram 


Function 


ERRCNT 

RDDST 

RDMAP 

WTRMAP 

PRERR 


Will ccmipute count error. 

Will read DSTTAP. 

Will read MAPTAP. 

Will write MAPTAP. 

Will print error computation. 
Will print classification map. 


PRMAP 


